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CHANNEL SIMULATION FOR DIRECT DETECTION 
OPTICAL COMMUNICATION SYSTEMS 

M. Tycz and M. W. Fitzmaurice 
Goddard Space Flight Center 


INTRODUCTION 

The rapidly increasing need for wideband communication systems for space-to-space and 
space-to-ground links in the late 1970’s and early 1980’s has stimulated an increased interest 
in optical communication technology. There are active NASA and DOD supported pro- 
grams which are directed toward developing high data rate, nominally 300 MBps to 1 GBps, 
optical communication systems for spaceflight use. A typical system might employ a 
Nd; YAG (neodymium doped, yttrium aluminum garnet) mode-locked laser transmitter 
using on-off-key (OOK), binary polarization (BPM) or quaternary modulation formats and 
a photoemissive direct detection receiver incorporating either high speed, gated photomul- 
tiplier tubes or avalanche photodiode/low noise, wideband amplifier combinations. Such 
systems are governed by Poisson statistics [ 1-2] . Their performance for OOK or BPM 
modulation formats and deterministic average receiver signal intensity have previously been 
analyzed and reported in the literature [3-4] . 

Optical communication systems look most attractive in the space-to-space link. In this con- 
figuration, it is easiest to take full advantage of their relatively small size and weight, their 
rugged and long-lifetime components and their low prime power requirements. However, it 
is in this configuration, for example low-earth orbiter to synchronous satellite link, that the 
operational system may require sub-arcsecond pointing of the laser transmitter beam. As a 
result, a dominant noise source may be the interaction of the transmitter far field irradiance 
profile and the instantaneous transmitter pointing direction (see Figure 1 ). Thus, the average 
receiver signal level may be random and obey beta statistics [5-7] . For the space-to-ground 
link the transmitter low power, light weight, etc. system parameters ean be taken advantage 
of, but the perturbing effects of the atmosphere must be realized (see Figure 1 ). Perform- 
ance analysis reported in the literature for the space-to-ground link incorporates the 
effects of atmospheric scintillation assuming a log-normal channel [8] . Although system 
performance analyses have not been completed for the Ricean and Rayleigh channels, 
theory does exist asserting that certain communication links will be subject to atmospheric 
scintillation resulting in time varying average signal amplitudes described by Ricean or 
Rayleigh probability distributions [9] . 

Coupling the advancement of the theoretical analyses of direct detection optical communi- 
cation systems for varied channels with the rapidly progressing systems development, it now 
becomes necessary to design a laboratory device for controlled testing of the communication 
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Figure 1. Direct Detection Optical Communication Channels 


system. It is required that the unit simulate the channel effects of the different link con- 
figurations. This document addresses the design and performance of one such Channel 
Simulator. 


CHANNEL SIMULATOR SYSTEM CONCEPT 


The purpose of the Channel Simulator is to provide the controlled effects of atmospheric 
scintillation and transmitter pointing inaccuracy for testing and evaluating direct detection 
optical communication systems for the space-to-space and space-to-ground links. As seen 
in Figure 2, the Channel Simulator transforms a constant signal 1 into a preselected random 
varying signal Ij^ . The probability distribution fi^ (i^ ) governing the modulated signal level 
of the Channel Simulator output 1^ (see Figure 3) will depend on the link to be simulated. 
The Channel Simulator described in this publication is designed to produce random signals 
whose time varying irradiance may obey log-normal, chi-squared, or beta statistics and whose 
time varying amplitude may obey log-normal, Ricean, or Rayleigh statistics. It consists of 
two major components: the linearized optical modulator and the processing electronics. Its 
transfer function is given by = C V where 1^ is the irradiance of the modulated light beam, 
C is a constant determined by the modulator crystal parameters and V is the voltage applied 
to the modulator. In general, V may be deterministic or random. If V is random with prob- 
ability density /y(u), then the probability density governing [10] will be 
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OPTICAL SUBSYSTEM 

The optical subsystem consists either of a linearized acousto-optic modulator or a linearized 
electro-optic modulator. One such system utilized by the experimenters was a KDP electro- 
optic modulator. Employed in the transverse mode of operation, the modulation field was 
applied normal to the direction of light propagation. Using the axes defined in Figure 4, the 
input polarizer was in the X'Z plane at an angle of 45° from the Z-axis. The output polarizer 
was rotated 90° relative to the input. The transfer function for this modulator system [11] 
is 
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where co is the frequency of incident light, c is the speed of light, and n^ are the ordinary 
and extraordinary indices of refraction, is the applied field, r^^ is the electro-optic coeffi- 
cient of the KDP crystal, and (n^ - n^) L^, <o/2c, the retardation due to natural birefringence, 
is a function of temperature and can be compensated. Therefore, consider 
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If we define = V/D and let 
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where is the voltage required for maximum transmission through the crystal. 
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If a series expansion of the right-hand side of Equation 5 is performed, then 
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A linearized optical modulator system can now be formed by simply inserting a square-root 
module between the applied voltage V and the modulator crystal, assuming operation in the 
range V < . 

The linearized acousto-optic modulator is designed using the same principal. Modulation of 
light due to the interaction of a light wave with a sound wave induces a diffracted light wave 
at the sum or difference frequency. The transfer function for a beam incident into a medium 
of a distance JZ [12] is given by 
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where is the wavelength of the incident wave in a vacuum, p is the photoelastic constant 
of the medium, Iacoustic acoustic intensity, is the velocity of sound in the acoustic 
medium, and p is the mass density. Defining Iq as a reference acoustic intensity yielding 
maximum deflection. 
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and noting that>/T is proportional to V, the applied voltage, then 
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Using the same expansion as before and assuming voltages V , then the transfer function 
becomes linear where once again a square-rooting device is inserted before the modulator 
crystal; 
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Results of this technique applied to the Channel Simulator will be presented later in this 
document. 


PROCESSING ELECTRONICS 

The analog processing electronics generate the random voltages necessary to drive the 
optical modulator for simulating the desired channels. As seen in Figure 5, it consists of 
Gaussian random signal generators fed into antilog converters, squaring modules, and 
similar devices which perform appropriate mathematical operations. Although certain 
operations are shared by more than one channel, the transfer function for each channel 
will be presented separately for clarity. 

Atmospheric Log-Normal Channel 

For the atmospheric scintillation channel log-normal mode, the random variable that obeys 
log-normal statistics is 1/1^, or the instantaneous irradiance divided by the mean irradiance. 
To simulate the log-normal behavior of this channel the following random variable needs 
to be generated: 
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where 2L is Gaussian with mean K and variance [13], The transfer function of the 
channel simulator log-normal channel (see Figure 6a) from the Gaussian input, monitored at 
test point TP7, to the modulated irradiance output is given by 
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where A is the log-normal offset, b is the gain established by the last amplifier stage, 

is a zero mean Gaussian random voltage with variance , and C is the linearized modu- 

A T p 7 

lator subsystem constant of proportionality. 

Redefining variables-in Equation 13 results in Equation 12, where 1^ = bC/10 and 
L=-(X.pp^ -I- A) InlO/IO. Additionally, when simulating the log-normal channel, the log 
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Figure 5, Channel Simulator Processing Electronics, Composite of All Channels 
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Figure 6. Channel Simulator Transfer Functions 































amplitude must have a mean equal to the negative of its variance [ 14] . Therefore, letting 
E |l} = - Var | L } = ~o^ , relationships are established for the values of A and c^xjp 7 


A = o^lO/ln 10 


(14) 


4 =a2 [10/ln 10]2. (15) 

The value of Ug is selected from atmospheric scintillation data and ranges nominally from 
o^ = 0.01 too2 = 0.1 [15]. 

Atmospheric Rayleigh and Ricean Channel 

For the atmospheric scintillation channel Rayleigh or Ricean mode, the random variable that 
is described by these statistics is the light amplitude. Until now, however, the linear transfer 
function of the Channel Simulator has been referenced to the irradiance level of the output. 
The light amplitude transfer function is thus 


A = 


(16) 


where C, = C^. 

A Rayleigh or Ricean random variable P can be generated by the following; 

P=[x2 + x|]’'^ (17) 

where for the Ricean case and Xj are independent Gaussian random variables with means 
Pj and pj variances Oj = a|. For the Rayleigh case, = Mj “ ■ 

The total transfer function of the Ricean* channel from the Gaussian inputs to the modu- 
lated light amphtude output is given by (see Figure 6b) 


P = 
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where Xpp 2 and X^-pj are zero mean independent Gaussian random variables with variance 
and d is the Ricean offset. Note also that for Rayleigh simulation, V and therefore the 
modulated irradiance is a chi-squared random variable with two degrees of freedom. 


*Since it is more general, we present this channel with the understanding that for the Rayleigh case the off- 
set is simply set to zero. 
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Free Space Channel 

For satellite-tosatellite optical communications, signal fades caused by pointing jitter may 
limit receiver performance. In this case, the modulation on the average received signal level 
will be described by beta statistics [5-7] . In simulating this fading condition the following 
must be satisfied: 


*0 ® 




(19) 


where 1^ corresponds to the peak of the spatial Gaussian irradiance profile, ^ is the transmitter 
beam half angle at the 1 /e^ points, and p is the beam displacement from the center of the 
receiving aperture. For typical systems, p^ is a chi-squared random variable of two degrees 
of freedom corresponding to two independent and orthogonal Gaussian tracking error 
signals. 

The total transfer function of the Channel Simulator in the beta mode of operation can be 
seen in Figure 6c. From the Gaussian random voltage inputs to the modulated irradiance 
output the relationship is 


I„ = lOCe-O" io/5)|(x2p2 + X2pj)/20I 
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where and are zero mean, equal variance, Gaussian random variables. Redefin- 
ing 

Iq = IOC 
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and 
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then Equation 19 is satisfied and the density of becomes 


.2 / , 
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The ratio of |/a is determined by both system design and the link to be considered and will 
vary typically from 3 to 10 [7] . 


COMPOWENT TRAiSiSFER FuTvCTiON TtSi KtsULTS 

Tests have been performed on the Channel Simulator processing electronics and optical 
modulator subsystem to verify the accuracy of the mathematical operations, the linearity of 
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the optical modulator, and the independence of the Gaussian sources. Figure 7 shows the re- 
sults of tests performed on the electronic modules which perform the required mathematical 
operation within the processing electronics. The equations representing the resultant curves 
were derived from the data assuming minimum mean square error fit to a straight line, least 
square fit to a power curve, and least squares fit to an exponential curve where appropriate. 

The transfer function of the linearized acousto-optic modulator system used by the authors 
was tested using a helium-neon laser input beam and an S-20 photomultiplier detector. Fig- 
ure 8 shows typical test results when the modulator bias and the photomultiplier tube voltage 
are set for linear operation. Since the slope of this transfer function is directly affected by 
the photomultiplier tube gain as well as the modulation index, this test was performed be- 
fore each of the channel simulations. 

As previously stated, operation of the Channel Simulator for the beta, Rayleigh, and Ricean 
modes assumes that the Gaussian random voltage inputs at test points TP2 and TPS are zero- 
mean, equal variance, and independent. Offset adjustments as well as variable gain stages are 
provided to ensure that the signals are zero mean and equal variance. Tests for independence 
were performed using the Hewlett Packard Correlator, Model 3729A. Recall that two random 
signals X and Y are independent if their cross-correlation coefficient is zero, that is. 


e{xy} -e{x} e{y} (23) 

VVar{x) Var{Y} 

where for this case, the means E{x} and e{y} are zero and the Var |X^ = Var {y| is 
simply the autocorrelation of X or Y evaluated at a time delay of zero, and E{xY} is the 
cross-correlation evaluated at a time delay of zero. The auto- and cross-correlations of 
Xjp 2 ^5 f^icean offset d = 0. Using these values at delay zero and 

substituting into Equation 23 resulted with the cross-correlation coefficient p = 0.006. 

STATISTICAL FUNCTION TEST RESULTS 

The usefulness of the Channel Simulator is measured in terms of its ability to accurately 
simulate the desired channel. Figure 9 illustrates the. performance of the log-normal 
channel simulation as compared to experimentally measured channel perturbations. The 
plot of log amplitude versus its cumulative probability is on a log-normal cumulative proba- 
bility graph. The solid line represents typical results from a GEOS-B satellite experi- 
ment [15] for which an argon laser beacon was transmitted through the atmosphere to a 
low earth-orbiting satellite. The instantaneous signal was detected onboard the spacecraft 
by a photomultiplier tube and was telemetered to ground for statistical analyses. The 
dashed lines represent data taken from the Channel Simulator system in which the signal 
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Figure 7 (continued). Transfer Functions of Mathematical Operations 



Figure 8. Transfer Function of the Linearized Acousto-optic Modulator 
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PERCENTAGE OF TIME THAT LOG - AMPLITUDE < X 
Figure 9- Comparison of Channel Simulator Log-normal Statistics with GEOS-B Data 


level was log-normally modulated. The linearity of the data indicates the modulated signal is 
indeed log-normally distributed and the slope specifies the strength of atmospheric scintilla- 
tion. The latter can be made to vary from a log-amplitude variance of 0.01 to 0.08. 

Performance of the beta channel requires operation of all the mathematical functions in the 
processing electronics. The output signal of the Rayleigh and Ricean channel is the input 
signal provided to the antilog converter for generating a beta-distributed random variable. 

Thus accurate simulation of the beta distribution necessitates good performance of the Ricean 
and Rayleigh atmospheric channel modes. 

Recall that the probability density function governing the beta mode of operation was given 
by Equation 22, Therefore, the cumulative density will satisfy 


r ^ .,- 7 . 0 , / 1 


14 




Taking the logarithm of both sides gives 


log (1m) ~^log lo- 

Plotting the logarithm of the cumulative beta density versus the logarithm of the relative 
modulated signal level irradiance results in a linear curve, the slope of which is the parameter 
of the beta distribution to be simulated. Figure 10 shows a typical plot of simulated beta 
modulation for which the ratio ^/a = 5.3. The Channel Simulator can be made to simulate 
the space-to-space link with beta perturbations ranging nominally from 1 < ^/a < 10. 

Although not directly applicable to simulating optical communication channel effects, an 
additional mode of operation has been coupled to the system. By adding an optional square^ 
rooting stage to the output of the Rayleigh and Ricean channel, the voltage output and, 
therefore, the modulated irradiance (as opposed to the modulated E-field) becomes Rayleigh- 
or Ricean-distributed. A test was performed on the operation of this Rayleigh random voltage 
output. Recall that a Rayleigh- distributed random variable P has cumulative density 

Fp(p)=-^/* (26) 

where is the parameter of the Rayleigh distribution. Performing the integration in Equa- 
tion 26 results with 


V-ln[l-F0t>)] 



(27) 


Figure 1 1 is a plot of Equation 27 using the square-root output voltage for data. The linearity 
of the data indicates that the voltage is Rayleigh-distributed and the slope indicates the value 
of the Rayleigh parameter, = 1 0.3 

CONCLUSIONS 

The system discussed in the previous sections for simulating the effects of space-to-space and 
space-to-ground optical communication channels is relatively easy to construct and enables 
controlled tests to be performed in the laboratory. The ultimate utility of the simulator is 
dependent on the validity of the theoretical models chosen to represent each channel. When 
coupled with a laser source, the Channel Simulator is capable of providing log-normal, Ray- 
leigh, Ricean, chi-squared, or beta irradiance modulation and voltages, as well as Rayleigh- and 
Ricean-modulated electric fields. It is intended that the simulator discussed in this document 
will be incorporated into a high datarate optical communication system for laboratory testing of 


15 






the composite system performance. Comparisons with computer-generated analyses will then 
be made to evaluate the accuracy of current analytical models for laser communication links 
undergoing fading. 
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